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Abstract 

A formulation is presented for modeling a resistive card in the 
context of the finite element method. The appropriate variational 
function is derived and for validation purposes results are presented 
for the scattering by a metal-backed cavity loaded with a resistive 
card. 
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1 Introduction 


A resistive card is an infinitesimally thin sheet of material which allows partial 
penetration of the electromagnetic field. Thin dielectric layers and very thin 
conductors whose thickness is less than the skin depth are examples of materials 
which can be modelled by resistive cards or sheets. Resistive cards are often 
used for radar cross section and RF power penetration control and as a result 
they have been studied extensively. Such studies have generally been done in 
the context of high frequency [1], [2] and integral equation solutions [3], [4], but 
to date the treatment of resistive cards within the context of the finite element 
method (FEM) has not been considered. Over the past few years, FEM has 
been applied to a variety of electromagnetic problems and it is thus important 
to incorporate the modeling of resistive cards in the FEM. In this paper we 
propose an FEM formulation which accounts for the presence of resistive sheets. 
To validate this formulation, results based on a physical modeling of the resistive 
sheet are also presented. In this case, the resistive sheet is equivalently replaced 
by a thin dielectric layer. The modeling of such a layer in the usual manner 
leads to larger and consequently inefficient linear systems, which is the primary 
reason for resorting to a mathematical modeling of the resistive sheet. Results 
based on the mathematical and physical modeling are presented in connection 
with the scattering by a metal-backed cavity in a ground plane and these are 
used to validate the proposed mathematical model. 

2 Formulation 

A resistive sheet is characterized by its resistivity R which is measured in Ohms 
per square. Mathematically, it satisfies the boundary condition [5] 

h T x (n r x E) = —R h r x (h + - H~) (1) 

where H* denotes the magnetic field above and below the sheet, E is the electric 
field and its tangential component is continuous across the sheet, and h r denotes 
the unit vector normal to the sheet pointing in the upward direction (+ side). To 
a first order, this boundary condition can be used to simulate the presence of a 
thin dielectric layer by setting [4], [5] 

*=]i^TTf (2) 
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In this, t is the thickness of the layer, Z 0 and ko denote the free space intrinsic 
impedance and wave number, respectively, and t T is the relative dielectric con- 
stant of the layer. Alternatively, a resistive sheet may be equivalently replaced 
by a thin dielectric layer having thickness t and a relative permittivity of 

~ = 1 - jZ ° 

6r Ic 0 tR' 

Generally, the accuracy of this simulation increases as the thickness t is decreased. 
Typically, t should not exceed one-tenth of the wavelength in the material. 

Let us now consider a finite element solution of the fields within a volume V 
subject to a given excitation. The volume consists of some inhomogeneous di- 
electric having relative permittivity and permeability e r and /r r , respectively, and 
we shall also assume that resistive cards may be embedded within the dielectric 
(see Figure 1). Within this volume, the electric field satisfies the vector wave 
equation 

V x ^-j-V xe|- k*e r E = -jk 0 Z 0 3 ial + V x (4) 

where (J int , M mt ) denote the impressed sources internal to the volume. Assuming 
a known tangential magnetic field on the boundary surface S 0 enclosing V, the 
electric and magnetic fields must satisfy the relation 

h 0 x V x E = —jk 0 ZoHrh 0 x H (5) 

in which n 0 is the outward unit vector normal to S a . This is simply a statement 
of Maxwell’s equations for the tangential boundary fields. 

For a finite element solution of the subject field, we introduce the functional 


( 6 ) 


F(E) "llllv [^ (V x E ) ' ( V >< E ) - • E dV 

+ JJJ v E • jk 0 Z 0 3 int - V x Q-M int j dV 

+^k Q Z 0 JJ s ^(n r x E) • (h r x E) dS 
+jkoZ 0 $ s E • (H x h Q ) dS 
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in which S r denotes the surface of the resistive sheet within V. We assert that 
F is stationary with respect to the true solution which satisfies (1), (4) and (5) 
and to prove this, we take the first variation of (6) with respect to E giving 


6F=JJJ v i(V x E) • (V x «E) - k\t r E ■ <5E 

+fv 6E -h 


jk 0 Z 0 3' nl — V x — M 


dV 

dV 


+jk 0 Z 0 JJ s ~(n T x E) • (n r x SE) dS 
-\-jkoZo (jj)$ £E • (H x h 0 ) dS. 


(?) 


On invoking a common vector identity, this can be rewritten as 

SF = jjj v jv x xeJ- kle r E + jk 0 Z 0 3 int - V x 


+ /l v [* Ex fe VxE ). 


dV 


+jk 0 Z 0 JJs j^(h r X E) • (h r x <5E) dS 
+jk 0 Z 0 (jJ s <5E • (H x h 0 ) dS 


■ 8E dV 


( 8 ) 


where we can now employ the divergence theorem for the second integral. Doing 
so and noting that the tangential magnetic field on S T is discontinuous we obtain 


1 


1 


6F= /// V x —V x E - kfc r E + jk 0 Z 0 3 mt - V x — M 


rint 


V 

+ 


•<5E dV 


IL 


s 0 +s++sr 


6E x 


1 


V x E 


• hdS 


+jk 0 Z 0 JJs -^(n r x E) • (h r x <5E) dS 
+jkoZo JJ <5E • (H x h 0 ) dS 


(9) 


where S + denotes the top (or + side) surface of the resistive sheet and S~ denotes 
its bottom surface. Noting that h r x E is continuous across S T and that the 
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normals on Sf and S T are in opposite directions, it follows that 


SF= Uf v | V x ( -j-V x EJ - klt T E + jk Q Z 0 J int - V x 


int 


•6E dV 


-jk 0 Z 0 JJ s [in, x (h T x E) + h r x (h + - H") • SE dS 

— H n 0 x { —V x E ] + jk 0 Z 0 h o x H 
JJSo L \fir J 


■ 6E dS. 


( 10 ) 


On invoking the conditions (1), (4) and (5), it is seen that SF = 0 proving our 
original premise that F is stationary with respect to the maxwellian field E. 

Having obtained the explicit form of the functional F , we may proceed to 
discretize it numerically in accordance with the finite element method by first 
subdividing the volume V into smaller volume elements and subsequently ex- 
panding the element field using vector basis functions. Substituting the field ex- 
pansion into the functional and applying the Rayleigh-Ritz procedure (equivalent 
to setting the first variation of F to zero), we can obtain a system of equations 
involving the interior and boundary electric fields and the boundary magnetic 
fields. For a unique solution of this system we must, however, specify a relation 
between the tangential electric and magnetic field which appear in the surface 
integral over S 0 . This can be provided by a boundary integral equation and the 
less accurate absorbing boundary conditions. If we assume the subject volume 
to be that occupied by the metal-backed cavity recessed in a ground plane, as 
shown in Figure 2, then S 0 reduces to the aperture area of that cavity. By in- 
voking image theory, the magnetic field on the aperture can then be expressed 
as 


H = H inc + H ref 


2jk 0 Y 0 JJ s 



Go (r, r') 


• [E(r) x z\ dS' 


( 11 ) 


where S a denotes the aperture surface, Go(r,r') is the free space scalar Green’s 
function, r specifies the observation point located on S a and I = xx + yy -f zz is 
the unit dyad. Also, H‘ nc denotes the magnetic field generated by sources in free 
space and H ref is the corresponding reflected field when the cavity’s aperture is 
shorted. A Galerkin’s discretization of (11) yields a system relating the aperture 
electric and magnetic fields and when it is combined with the finite element sys- 
tem we can obtain a complete system for a unique solution of the interior and 
boundary fields. The system obtained in such a manner will be partly sparse 
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and partly full. In particular the volume integrals and that over S T in (7) lead 
to a sparse submatrix involving the interior fields of the cavity. However, the 
discretization of (11) over S a renders a full Toeplitz submatrix since the bound- 
ary integral is convolutional. Consequently, by resorting to an iterative solution 
such as the conjugate or biconjugate gradient method in conjunction with the 
fast Fourier transform, the need to generate the Toeplitz matrix is eliminated 
thus maintaining the O(N) storage requirement, characteristic of finite element 
solutions. The details pertaining to this implementation are discussed in [6]- [8]. 
In the next section we only present some results aimed at evaluating the accuracy 
of the proposed resistive sheet model. 


3 Numerical Results 

Let us consider the metal-backed rectangular cavity illustrated in Figure 2. The 
cavity is assumed to be empty (no internal sources) and is illuminated by a plane 
wave in the <j) = 0 plane. To implement the aforementioned solution, the cavity 
is subdivided into rectangular bricks and the results of the solution are shown 
in Figure 3. These are radar cross section (RCS) patterns and refer to a 1A 
deep cavity whose aperture is also 1A x 1A. The RCS pattern in Figure 3(b) 
applies to the cavity which is loaded with a resistive sheet of 100f2/n placed at 
its aperture, whereas the result in Figure 3(a) is for the untreated empty cavity. 
The simulation of the resistive sheet was done through direct discretization of the 
first variation of the functional F as given in (6) and alternatively by modeling 
the resistive sheet as a dielectric layer of thickness A/20 having the dielectric 
constant computed from (3). As shown in Figure 3(b) the results based on the 
two simulations are in reasonable agreement and the differences among them is 
due to the finite thickness which was necessarily introduced in the physical model 
of the resistive sheet. To support this statement, in Figure 4 we show the RCS 
computed at 6 = 40° for several different thicknesses of the dielectric layer. The 
RCS value corresponding to the ideal resistive sheet simulation coincides with 
t = 0 of the plot. As seen, for smaller t the physical simulation is in better 
agreement with the ideal resistive sheet simulation. When t = A/20, one can 
expect a difference of one to two dB between the physical and mathematical 
simulation. 

The second geometry which was considered is a circular metal-backed cavity 
again situated in a ground plane. In this case the aperture of the circular cavity 
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is loaded with a sheet having non-uniform resistivity given by 


Zo 


R(p) = 


o.i + iofV^ 


t CO 


p > 0.5a 
elsewhere on S a 


( 12 ) 


where a denotes the radius of the aperture. Results with and without resistive 
loading for a cavity having a = 1 in. and a depth of 0.25 in. are shown in 
Figure 5. These RCS patterns were computed at 16GHz and the incident field 
was a plane wave polarized along the 6 or <f> directions. Again, the data in Figure 
5 demonstrate the validity of the proposed mathematical model. Also, as in 
the case of the rectangular cavity the presence of the resistive cards reduces the 
overall RCS of the cavity and this reduction is primarily due to the reduced field 
intensity near its perimeter. 


4 Conclusions 

A variational formulation was presented for modeling resistive cards within the 
context of the finite element method. Results based on the discretization of the 
variational expression were also presented and these were aimed at demonstrat- 
ing the validity and accuracy of the proposed mathematical model. Modeling 
of conductive sheets [5] can be accomplished by employing the magnetic field 
formulation yielding a functional which is the dual of (6). However, when the 
conductive and resistive sheets are both present within the computational do- 
main one cannot avoid the introduction of double nodes on both sides of the 
surface occupied by the conductive or resistive sheets. If an electric field formu- 
lation is employed as given here, double nodes must be placed on the surface of 
the conductive sheet. The functional in (6) must also be supplemented with the 
additional integral 

J -k 0 Z 0 j] s R m [h c x (E + - E")] • [h e x (E + - E-)] dS (13) 

where S c is the surface occupied by the conductive sheet, Rm denotes its conduc- 
tivity and n c is the normal to the sheet. Other than the requirement to introduce 
double nodes on S c , the implementation of (13) is straightforward. 
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FIGURE CAPTIONS 


Fig. 1 Cross section of a dielectric volume enclosing a resistive card. 

Fig. 2 Geometry of a cavity-backed aperture in a ground plane. 

Fig. 3 Monostatic RCS for a square aperture (1A x 1A) backed by a rect- 
angular cavity (1A x 1A x 1A) in the <f> = 0 plane, (a) Empty cavity 
without resistive loading, (b) Empty cavity whose aperture is covered 
with a resistive sheet of a resistivity of 100 f !/□. (Solid and dashed 
lines correspond to results based on the mathematical simulation; cir- 
cles and squares refer to results based on the physical modeling of the 
resistive sheet). 

Fig. 4 Monostatic RCS at (0 = 40 °,(f> = 0) of the square aperture in Fig. 
3(b) as a function of the laryer thickness simulating the resistive sheet 
at the aperture. The e r of the layer is given by (3) with R = 100U/Q. 

Fig. 5 Monostatic RCS for a circular aperture (1 inch in diameter) backed 
by a circular cavity (1 inch in diameter and 0.25 inch in depth) at 
16 GHz. (a) Empty cavity without resistive loading, (b) Empty cav- 
ity covered with a resistive sheet whose resistivity is given by (12). 
(Solid and dashed lines correspond to results based on the mathemati- 
cal simulation; circles and squares refer to results based on the physical 
modeling of the resistive sheet). 
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Fig. 4 















